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ON THE LAGRANGIAN DESCRIPTION OF UNSTEADY 
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Leon L. Van Dommelen* and Stephen J. Cowleyt 
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Lewis Research Center 
Cleveland, Ohio 44135 

Summary 

Although unsteady, higli-Reynolds-number, laminar boundary layers have convention- 
ally been studied in terms of Eulerian coordinates, a Lagrangian approach may have sig- 
nificant analytical and computational advantages. In Lagrangian coordinates the classical 
boundary-layer equations decouple into a momentum equation for the motion parallel to 
the boundary, and a hyperbolic continuity equation (essentially a conserved Jacobian) for 
the motion normal to the boundary. The momentum equations, plus the energy equation if 
the flow is compressible, can be solved independently of the continuity equation. Unsteady 
separation occurs when the continuity equation becomes singular as a result of touching 
characteristics, the condition for which can be expressed in terms of the solution of the mo- 
mentum equations. The solutions to the momentum and energy equations remain regular. 
Asymptotic structures for a number of unsteady three-dimensional separating flows follow 
and depend on the symmetry properties of the flow (e.g. line symmetry, axial symmetry). 
In the absence of any symmetry, the singularity structure just prior to separation is found 
to be quasi two-dimensional with a displacement thickness in the form of a crescent shaped 
ridge. Physically the singularities can be understood in terms of the behavior of a fluid 
element inside the boundary layer which contracts in a direction parallel to the boundary 
and expands normal to it, thus forcing the fluid above it to be ejected from the boundary 
layer. 
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1. Introduction 


A major feature of unsteady large- Reynolds-number flow past a rigid body is the shed- 
ding of vortices from the surface of the body. Such vortices alter the forces exerted on the 
body dramatically (McCroskey & Pucci 1982). A more complete theoretical understanding 
of vortex shedding would be advantageous in the design of air, land and water transport. 
Theoretical models of vortex shedding also have application, inter alia , in the description 
of air flow over hills and water waves, water flow over sand ripples, and blood flow through 
curved and constricted arteries and veins. 

A classical example of vortex shedding develops when a circular cylinder is set into 
motion in the direction normal to its axis. This example was first studied by Prandtl (1904), 
and the process by which an initially attached boundary layer develops into a separated flow 

with detached free shear layers has been clearly illustrated by the experiments of Nagata, 
Minami & Murnt.a (1970), and Bouard &■ Ooutanreau (1980). The term ‘Reparation’ will 
in this paper be used to refer to the ‘breakaway’ of a thin layer of vorticity from the 
surface of a body. This definition of separation is close to that of both Prandtl (1904) 
and Sears & Telionis (1975). In particular, Sears &: Telionis speak only of separation when 
the penetration of the boundary-layer vorticity away from the wall becomes too large to 
be described on the usual 0(Re ~ *) boundary -layer scale (Re is the Reynolds number of 
the flow, and is assumed large). Therefore, once separation has developed the classical 
attached flow solution will, in general, no longer be valid. 

The first theoretical advance in understanding the unsteady cylinder flow at high 
Reynolds numbers was made by Blasius (1908). He explained the occurrence of flow 
reversal inside the attached unsteady boundary layer which is set up immediately the 
cylinder starts to move. In the case of steady flow past a rigid surface, flow reversal is 
often accompanied by separation. However, Moore (1958), Rott (1956) and Sears (1956) 
all realized that zero wall shear is not necessarily related to separation in unsteady flow. 
Sears & Telionis (1975) noted subsequently that their definition of separation is consistent 
with the termination of the boundary-layer solution in a singularity. Such a singularity 
will be referred to as the separation singularity, and the time at which it develops as the 
separation time. 


2 



A considerable number of numerical computations have attempted to verify the ex- 
istence of a singularity in the boundary-layer solution for the circular cylinder problem. 
The first convincing evidence that a singularity forms within a finite time was given by 
Van Dommelen & Shen (1977, 1980a, 1982). In a Lagrangian computation, with fluid par- 
ticles as independent coordinates, they found that a separation singularity develops after 
the cylinder has moved approximately | of a diameter. The existence of this singularity 
has been confirmed by the finite difference numerical calculations of Ingham (1984) and 
Cebeci (1982) (however see Cebeci, 1986), and the computer extended series solution of 
Cowley (1983). These calculations were all based on Eulerian formulations. A similar 
two-dimensional separation singularity has been observed using Lagrangian procedures on 
an impulsively started ellipse at several angles of attack (Van Dommelen, Wu, Chen k 
Shen, unpublished results), on airfoils (Wu 1988), in turbulence production (Walker 1988), 
on an impulsively started sphere (Van Dommelen 1987), and using an Eulerian scheme in 
leading edge stall (Cebeci, Khattab k Schirnke, 1983). 

Excluding vortex methods, flows with free surfaces, and some more specialized com- 
pressible flow computations, Lagrangian coordinates have not been as widely used as their 
Eulerian counterparts in fluid mechanics, especially for boundary-layer flows. Yet for some 
flows, such as unsteady flows in which advection dominates diffusion, Lagrangian coordi- 
nates seem more appropriate (e.g. see the inviscid calculations of Stern k Paldor (1983), 
Russell k Landahl (1984) and Stuart (1987)). As far as unsteady separation is concerned, 
the advantage of a Lagrangian approach stems from the fact that in these coordinates 
the classical boundary-layer equations decouple into a momentum equation for the motion 
parallel to the boundary, and a continuity equation for the motion normal to the boundary 
(Shen 1978). The solution of the former equation can be found independently of the latter. 
Moreover, while the time that the separation singularity develops can be identified from 
the solution to the momentum equation, only the solution to the continuity equation is 
singular (see section 2). 

An important consequence of the Lagrangian approach is that simple descriptions can 
be found to a wide variety of separations in one-, two- and three-dimensional unsteady 
flows. In this paper we consider unsteady flows in general, then in part 2, (Van Dommelen 
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1989), the separation process that occurs at the equatorial plane of a sphere which is set 
into a spinning motion is examined in detail. 

In the next section we develop the simple analytic machinery needed to find self- 
consistent three-dimensional separation structures for both compressible and incompress- 
ible fluids. Some of the properties of the Lagrangian version of the boundary-layer equa- 
tions are also discussed. In section 3 the Lagrangian structure for three-dimensional sep- 
aration is derived under the assumption that the flow can be completely general, then in 
section 4 the changes in structure are discussed when various symmetries restrict the flow 
geometry. 


2. Lagrangian Formulation 

The Lagrangian description of boundary-layer flow uses fluid particles (i.e. infinitesi- 
mal masses of fluid) as the basis of the coordinate system. A convenient coordinate system 

for the fluid particles (£,T7,C) is given by the initial Eulerian position of the particles (see 
Lamb (1945) for example): 

£ = (£, * 7,0 = (*,!/> z ) at < = 0 • ( 2 . 1 ) 

The precise form of the Lagrangian solution depends on the particular reference time, 
defined here as the start of the motion, but the physical solution is independent of it. 

Following Rosenhead (1963) we assume that the position coordinates * and z describe 
an orthogonal coordinate system on the surface of the body in question. The lengths of the 
line elements dx and dz are taken as hi dx and h 3 dz respectively. The coordinate normal 
to the surface is denoted by y, which is scaled with the square root of the reference shear 
viscosity. 

In Lagrangian coordinates, conservation of volume for a compressible fluid can be 
expressed in terms of a Jacobian determinant as follows (e.g. Hudson 1980): 

pH{x,z)J(x,y,z) = poH 0 , (2.2a) 
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where 


J(x,y,z) 


x 





x ,< 

y,v 

y,< 

z ,v 

z ,< 


» po{Z,v,0 = p(t,v,(,o) , 


(2.2b, c) 


H(x,z) = hi(x,z)h 3 (x,z) , H 0 = H(£, <) , (2.2d,e) 

* s the density of the fluid, and a subscript comma denotes a Lagrangian deriva- 
tive. The velocity components of the flow are related to the fluxions of position by 


u = hi(x,z)x , w = h 3 (x,z)z , 


(2.3a, b) 


where a dot represents a Lagrangian time derivative. 

For compressible flow, the momentum and energy equations are (e.g. Rosenhead 
1963): 

... w 1 

(2.3c) 


w 1 

p(u + ( uh lz - wh 3x ) — ) = p x + D y (nD y u) + pg x , 


'll \ 

p(w + ( wh 3x - uh iz ) — ) = p z + Dy(pDyw) + pg z , 
P^P + 9 Tp P ~~p 9 ^ P ^ DyU ^ + ( D v w ) 2 ) + D y(^ D y T ) , 


(2.3d) 


(2.3e) 


where p is the scaled shear viscosity, a is the Prandtl number, and g x and g z are the 
components of the acceleration of gravit}'. The temperature, T, and internal energy, e, are 
assumed to be functions of density and pressure, while the pressure, p, is a knowm function 
of x, z and t; thus 

u w 

P = pt + J-Px + T-Pz • (2-3/) 

rl\ rl 3 

For an incompressible flow p = 0 and e is taken to be a function of T and p. 

Although the y-derivative D y is Eulerian in nature, it can be written in the Lagrangian 
form (see also Shen 1978): 

p(t,V, (,t)H{x,z)J{x,u 1 z) 


Dy U = 


From (2.4) it follows that at a fixed wall the Eulerian D y and Lagrangian d/ dy operators 
differ only by the density ratio, which leads to simplifications in the calculation of the wall 
shear. 
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Allowing for a moving boundary, appropriate boundary conditions to (2.3) are: 


(u,w,p) = (u b (x,z,t),w b (x,z,t),p b (x,z,t)) on y = 0 , (2.5a) 


(u,w,p) -*• (u e (x,z,t),w e (x,z,t),p e (x,z,t)) as y ->• oo , (2.56) 


where u b and w b specify the velocity of the boundary in the x- and 2 -directions respectively, 
u c and w e are the corresponding external slip- velocities, p e is the external flow density, 
and the wall density p b can be given implicitly as the temperature at the wall. Ordinarily, 
these boundary conditions translate immediately to the Lagrangian domain by means of 
(2.3a, b). In the case of suction or blowing through the wall, they must be applied at an q- 
boundary moving through the Lagrangian domain, however, the wall boundary conditions 
turn out to be of little importance for the local analysis of this paper. 


The principle advantages of Lagrangian coordinates derive from the absence of both 
the normal particle position y and the normal velocity component v from (2.3) and (2.4). 
Consequently, the particles’ motion, as projected onto the surface of the body (x, 2 ), can 
be found independently of the normal particle position y. Subsequent integration of the 
Jacobian (2.2) along lines of particles at constant projected position (x, z) yields the normal 
particle position 

f po Ho ds (n a\ 

.mr 7 ~ a fT-i ’ ( 2 ‘ 6 ) 

Jo 


pH |Vx A V«| 


where ds 2 = d£ 2 + drj 2 4- d£ 2 , Vx = (x t ( ,x >v ,x t () is the Lagrangian gradient, and the 
integral is performed in the Lagrangian £;t) coordinate system along the lines of 

constant x, z, and t, i.e. lines which in physical space are vertical through the boundary 
layer. 

The central issue of this paper can now be stated: we hypothesize that during the evo- 
lution toward separation, the projected position (x, z) can remain regular, and commonly 
does remain regular. When true, such regularity strongly restricts the possible behavior of 
x and z near separation, and to characterize separation we need only identify the nature 
of solutions to the continuity equation (2.2) or (2.6) - an equation which is much simpler 
than the momentum equations. The remaining ambiguity in the behavior of x and z is 
resolved using arguments of symmetry. 
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Various arguments to justify our hypothesis can be given. One of them is self- 
consistency. If it is assumed that x, z,u,w, and p are non-singular at the separation time 
t 9 , then the solution to the Lagrangian momentum equations can be expanded in powers 
of (t — t s ) to any algebraic order. In contrast, the usual Eulerian asymptotic expansions 
show only that the first few terms in the expansions axe self-consistent. 

As another argument, Van Dommelen (1981) showed analytically that the inviscid, 
incompressible two-dimensional equations have solutions, x, z, which are regular func- 
tions of the Lagrangian variables, although y(£,t) is singular (this analysis can be further 
developed by expanding in powers of a small coefficient of viscosity). Yet this example is 
somewhat artificial; physically it would require that during the evolution of the boundary 
layer the coefficient of viscosity was changed significantly by some external means. 

A more powerful argument is possibly the capability of the analysis in this paper to 
reproduce and extend several known separation processes previously analyzed in Eulerian 
coordinates. However, the most convincing argument is provided by actual numerical solu- 
tions of the Lagrangian boundary-layer equations. For example, Van Dommelen & Shen’s 
(1980a) computation of the boundary layer on an impulsively started circular cylinder 
provided direct numerical evidence as to regularity of the momentum equation. Further, 
it is in remarkably close agreement with the results obtained by Cowley (1983) using a 
series extension technique. In particular, Cowley (1983) finds a singularity in the solution 
at the same time and position as the Lagrangian computations. Ingham (1984) performed 
an Eulerian Fourier series expansion of the solution in the direction along the cylinder. 
By carefully increasing the order of expansion as the spectrum expands due to the incipi- 
ent singularity, he obtained results in close agreement with those of both Van Dommelen 
& Shen (1980a) and Cowley (1983). The fact that these three very different procedures 
were found to produce results in excellent agreement with one another until very close to 
the breakdown of the solution at separation is reassuring, since a number of more con- 
ventional finite difference computations (e.g. Telionis & Tsahalis (1974), Wang (1979), 
Cebeci (1986)) give significantly different results. Yet the results of Henkes & Veldman 
(1987) remain in agreement with the three unconvential methods until relatively close to 
the singularity, but disagree with Cebeci (1986) at a significantly earlier time. One of the 
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difficulties with conventional finite difference procedures, as pointed out by Cebeci (1986), 
is the need to satisfy the CFL condition, a condition which is implicitly satisfied by the 
three procedures of Van Dommelen & Shen, Cowley and Ingham. 

Clearly in any numerical Lagrangian computations, it is not possible to prove that the 
solution is regular, since the inevitable upper limit on resolution means that high order 
singularities are difficult to resolve. However, in the accompanying numerical study, part 
2, the boundary layer at the equatorial plane of a spinning sphere is solved using up to 
1000 mesh points across the boundary layer. Even at such high resolution, no trace of 
singular behavior was observed, and derivatives of high order could be evaluated precisely. 

When the fact that solutions to the momentum equations are regular is accepted, 
(and for compressible flow in addition the density must be regular), the next question to 
arise is what implications such regularity has for the structure of the separation process. 
First, only the continuity equation can develop singular behavior, and from (2.2) or (2.6) 
it follows that this is only possible if the Lagrangian gradients of * and z become parallel, 
i.e. if at some point s 

Vx = A 4 V 2 , (2.7a) 

where A 4 is a constant. Generally, the point s of interest is the particle and time at which 
(2.7 a) is satisfied for the first time. The condition (2.7 a) is a three-dimensional extension 
of the two-dimensional condition first pointed out by Van Dommelen & Shen (1980a); it 
requires that a Lagrangian stationary point, Yn = 0, exists for an oblique coordinate 

n = x — X s z . (2.76) 


An alternate way to phrase the condition for singular y is to define a unit vector in 
n-direction tangential to the wall, 

n = (n s ,n z ) = %== , (2.7c) 

V 1 + 

in which case a singularity occurs when, for all infinitesimal changes d£ in fluid particle, 

n • dx = 0 , dx = (dx,dz) . (2.7 d,e) 
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This implies that an infinitesimal particle volume d£dr]d( around point s has been com- 
pressed to zero physical size in the n-direction. But since particle volume (or mass in 
compressible flow) is conserved, this compression in the n-direction along the wall is com- 
pensated for by a rapid expansion in the y-direction (see figure 1), which drives the fluid 
above the compressed region d(dr)d( ‘far’ from the wall to form a separating vorticity 
layer. 

From (2.6) it can be shown that this process constitutes separation in the sense of Sears 
k Telionis (1975), since the particle distance from the wall becomes too large, ‘infinite’, to 
be described on the usual boundary-layer scale. Note that the assumed regularity of x and 
z does not allow an infinite expansion in the direction parallel to the wall but normal to n; 
the particle can only expand strongly in the direction away from the wall. Similarly for a 
compressible fluid, the assumed regularity of p is inconsistent with an infinite compression 
of the particle volume. (At present there is no direct numerical evidence for the regularity 
assumption in the compressible case, although it is of course self-consistent). 

From (2.7) we can derive generalized so-called Moore- Rott- Sears (MRS) conditions at 
the stationary point, similar to the conditions formulated by Sears k Telionis (1975) for 
two-dimensional flow. The form of the Eulerian D y operator (2.4) implies using (2.2 b) and 
(2.7a) that the vorticity vanishes at that point, i.e. 

D y u = D y w = 0 at Vn = 0 . (2.8a) 

In fact, the D y -operator vanishes for all quantities which remain regular in the Lagrangian 
domain. 

The second Moore- Rott-Sears condition is more complicated. Since (2.7a) is equivalent 
to two conditions on the Lagrangian derivatives of x and z , in three dimensional space we 
expect it to be satisfied on a curve of particles for times beyond the first occurrence of 
separation (c.f. section 3 and sub-section lc). The Eulerian projection of the singular 
curve on the wall will be denoted by x.mrs = ( x MRS,z\fRs) and the Moore- Rott-Sears 
condition concerns the motion of this projected curve. To derive it, we focus attention on 
an arbitrary point s on the singular curve (rather than our usual choice in which s is the 
first point at which a singularity occurs). First we consider a Lagrangian differential 
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along the singular curve passing through point s, keeping time constant. Since x and z 
are functions of £ and t only, d£ corresponds to a change in Eulerian position along the 
projected curve which satisfies (2.7d), 

n • dx-MRS = 0 , (2.86) 


so that the singular curve is normal to the local vector n. As for any curve, the propagation 
velocity of this curve is given by the component of the propagation velocity of points on 
the curve in the direction normal to the curve. To find an expression for it, we now 
consider a total differential in Lagrangian space-time at the point s, resulting in changes 
dxMRS = dxMRS + *sdt and dzMRS = Ozmrs 4- z a dt. Since ( Oxmrs, Qzmrs ) satisfies 


(2.8b), 


n 


dxMRS 

di 


n • umrs 


umbs = { 6 », Zs ) 


(2.8c, d) 


which shows that the propagation velocity of the singular curve equals the flow velocity of 
the singular particle s at the considered position ( *mrs,zmrs )■ 

While this three-dimensional form of the MRS condition seems new, the general ap- 
plicability of the two-dimensional case is fairly well established both theoretically (Moore 
1958, Sears Sc Telionis 1975, Williams 1977, Shen 1978, Sychev 1979, 1980, Van Dorarae- 
len Sc Shen 1980b, 1982, 1983a, b, Van Dommelen 1981) and experimentally (Ludwig 1964, 
Didden Sc Ho 1985). 

We can also verify the notion of Sears Sc Telionis (1975) that unsteady separation 
occurs in the middle of the boundary layer rather than at the wall. In the absence of a 
transpiration velocity, the motion of points on the wall equals the motion of the boundary- 
layer particles at the wall, cf. (2.5) and (2.3a, 6). Thus a fluid particle at the wall can only 
contract to vanishing size in the n-direction if the wall itself performs the same contraction, 
which is not possible for a solid wall. 


In the next sections the nature of the separation process is analyzed. First we form 
local Taylor series expansions for the regular solutions to the momentum equations near 
the stationary point, and then we expand the solutions of the continuity equation in an 
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asymptotic series. This procedure is similar to the one followed by Van Dommelen & 
Shen (1982) for two-dimensional separation. In contrast to the steady viscous singularities 
of Goldstein (1948) and Brown (1965), and the ideas of Sears &: Telionis (1975), the 
unsteady singularity is essentially inviscid in character and consists of two vortex sheets 
separated by an increasingly large central inviscid region (as found by Ockendon (1972) 
for a rotating disc with suction, and by Sychev (1979, 1980), Van Dommelen & Shen 
(1980b, 1983a, b), Williams & Stewartson (1983) and Elliott, Cowley & Smith (1983) for 
steady separation over up- and down-stream moving walls). The leading order asymptotic 
structure of the unsteady singularity has also been recovered by Van Dommelen (1981) as 
a matched asymptotic solution to the Eulerian boundary-layer equations. More generally, 
Elliott et al. (1983) showed that there is a certain amount of arbitrariness in the Eulerian 
expansions. The Lagrangian expansion resolves such arbitrariness by the assumption, 
(supported by various numerical data, see Van Dommelen & Shen (1982), the closing 
remarks of subsection 4c, and part 2), that the leading order coefficients in the Taylor 
series expansion near the stationary point are non-zero. 

3. Three-dimensional separation singularities 

In this section we find the leading order term of an asymptotic analysis which describes 
the local structure of the flow as unsteady separation is approached. The time and position 
at which the separation singularity first develops will be denoted by the subscript s, thus 
for example 

(Vn), = 0 , (3.1a) 

where n is the oblique coordinate corresponding to the initial separation, defined in (2.76) 
as 

n = x — X s z . (3.16) 

Note that the definition of the x- and z-coordinates can simply be interchanged if n and 
z are not independent coordinates. In index notation, (3.1a) can be written as n* = 0, 
where we will adopt the convention to omit the subscripts comma (to indicate Lagrangian 
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derivatives) and s (to indicate the separation particle at the separation time) if they occur 
together (i.e. rii = (n ; ) g ). 

The solution of the continuity equation (2.2) for y can be greatly simplified by a 
number of coordinate transformations for both the particle position coordinates ( x,z ) and 
the Lagrangian coordinates (£,q, C)- Here we will select transformations which preserve the 
Jacobian J (2.2b), since these are algebraically more simple than transformations which 
preserve the physical volume HJ, or mass pHJ. 

As a first transformation, we drop the position coordinate x in favor of n, shift the 
Lagrangian coordinate system to the separation particle s, and rotate it, resulting in the 
set of coordinates 


n = x — \ s z , z , 


k ~ X] a *v(£j — * 


(3.2a, b, c) 


where a,j is an orthonormal rotation matrix which is is chosen to eliminate the mixed 
derivatives n J2 , n J3 , and n 2 3 - Therefore, expanding n and 2 in a Taylor series expansion 
about the separation point, we obtain 

3 / 3 


n — n, + ^ \nakl + . . . + St I n, + ^ n< fcj + . . . j + . . . , 

»=i ^ ' 


(3.3a) 


3 

z = z, + ^ Z{ki - f . . . + 6tz s , (3.36) 

t=l 

where St — t — i s . 

However, if t s is the first time that a stationary point occurs, the Taylor series co- 
efficients in (3.3) cannot be completely arbitrary: the singularity condition may not be 
satisfied anywhere for St < 0. The condition for a singularity to exist for earlier times at 
some neighboring point is, in terms of n and 2 , 


n,i - (A - A 4 )z )f ; = 0 


(3.4a) 


where A is the ratio between V 2 and Vx at the neighboring singular point. Expanding 
(3.4a) in a Taylor series, we obtain 

nuki — ZiSX + riiSt + . . . = 0 fort = 1,2,3 , (3.46) 
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where 6A = A — A,. If all three coefficients Tin, *i22> and ^33 were non-zero, a solution 
to (3.46) could be found for St < 0, contradicting our assumption that the singularity 
develops first at St = 0. (Strictly, because of the higher order terms omitted in (3.46), 
the solution must be found iteratively, however, the iterations converge because the higher 
order terms act as a contraction mapping sufficiently close to point s). Therefore at the 
first occurrence of separation, at least one of tin, 7122 ? or n 33 roust be zero, and we will 
reorder (ki,k 2 y k 3 ) such that nn vanishes. In addition, the coefficients Ti 22 ,ti 33 ,zi cannot 
all be non- zero, since by solving for S\, & 2 , and k 3 , it again follows that a singularity exists 
for St < 0. Without loss of generality, we assume that zj is zero, since if either n 2 2 or n 3 3 
vanishes, the (ki,k 2 ,k 3 ) coordinate system can be rotated further to eliminate z \ . 

It follows that in some suitably oriented Lagrangian coordinate system the conditions 
n n — Z\ =0 are necessary at the time when separation starts. This implies two addi- 
tional conditions on *(£,<) and z(£,<), besides the two conditions implicit in (3.1a). Since 
Lagrangian space-time is four-dimensional, in general we do not expect that more than 
four conditions can be satisfied at any time. Hence, in the remainder of this section we 
will assume that the values of the remaining derivatives can be completely arbitrary and 
in general non- zero. 

However, when the functions x and z are not arbitrary, but restricted by constraints of 
symmetry in the flow, the latter assumption needs to be reconsidered, since the symmetry 
requires that various derivatives must vanish. Examples are two-dimensional flow, and the 
flows discussed in the next section. 

Under the assumption that the remaining coefficients in the Taylor series have arbi- 
trary values, the transformation 


z — z ^(£ 4 ,/) , ft — 71 n(£ a ,t) A ^ ^ z , 

7 , 7 TI22 Z3 ~ U33 Z2 k 3 7 TI33 Z2 Al2 4 - TI22 23 ^3 

«1 = «J , *2 = 7===========- , fc 3 - - ==. = ==- 

V n 22 Z 3 d- n 33 2 2 


where 


\ n 22 2 3 + n 33 Z 


A< 2 > = 


7122^33 


3 2( 7122^1 + n 33 z\) 

eliminates the n 33 term. The final coordinate transform 


/] — k\ + — k 3 , I2 — k ?2 , l 3 — k 3 , 

n m 


(3.5a, 6) 
(3.3c, d,e) 

(3.5/) 


(3.6a, 6, c) 


13 




z = z 


(3.6d,e) 


n = n — A^z 3 — p a 6iz 


where 

\(3 ) ”iii ”333 — 3n in ^113^33 + 2n 3 13 n in n 3 — n u3 hi pict \ 

A » - ^2 ~3 » / i * — - - > 

on ni z 3 ^111 *3 

diminates the n H 3, n 333 , and n 3 derivatives. 

The transformed position coordinate n corresponds to an oblique coordinate measured 
from a moving, curved line through the separation particle, viz. 


where 


n = x — ®o(z,t) , (3.7a) 

i = , z = z-z(£,,t) , (3.76, c) 

x 0 (z,t) = A,z 4- A^z 2 4- A^z 3 + p,8tz . (3.7d) 


Note that the curved line x = xo(z,t), which can be viewed as the line along which 
the separation initially develops (see below), does not have a singular shape at the first 
occurrence of separation. 

The Taylor series expansions for n and z near the separation point become 


n - Iriaa/j + \n*jkkljlh + ••• + St'^nik + .. . , (n m = n 333 = n 3 = 0), (3.8a) 

ijk i 

z = Z 2 I 2 + 23/3 + .. . . (3.86) 

The diaracteristics of the Jacobian equation (2.2) for y are, in terms of the new coordinates, 

~ T ~ = ' P ^{-^3^22^2 + ••• } > ( 3 . 9 a) 

ay poTio 

~J~ = ~ 77 {-Z3 ( 3 TZH 1 Zj + |n,i 33 / 2 + TliSt) + . . .} , (3.96) 

dy p 0 H 0 

= — —{ — ^2(1^111^1 + 2 n i33^ 3 + TliSt) + .. .} , (3.9c) 

dy po H 0 

with a singularity occurring when all three right hand side expressions vanish (note that 
not all three are independent). Near the point s, (3.9a) is zero on a surface approximating 
the I 2 — 0 plane, while both (3.96) and (3.9c) vanish at points depending on the nature 
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of the quadratic expression (|nm^i + 5^133^3)- If this quadratic is hyperbolic, singular 
particles occur along hyperbolic lines regardless of the sign of St. Thus, if St = 0 is to be 
the first time that separation occurs, the quadratic must be elliptic, and of the same sign 
as the constant term when St < 0. This requires n m ni 3 3 > 0 and iJniiii < 0; we will 
choose the positive / x -direction such that 

n 22 n lll > 0 , n 22 n 133 >0 > n 22 n l < 0 • (3.10a, 6, c) 


The Lagrangian description of the separation process can now be completed by the 
determination of y at times shortly before the initial occurrence of separation. At t = t a 
the boundary-layer approximation is obviously no longer valid because from the integral 
(2.6) it follows that y becomes infinite at the stationary point. However, the rate of growth 
near this point can be found by means of an asymptotic expansion. To find local scalings, 
we follow the guiding principles of Van Dyke (1975). In general, we attempt to scale the 
Lagrangian coordinates li and the position coordinates n, z and y to variables Li,N,Z , 
and Y such that in the inner region the Jacobian equation for Y, i.e. 


Jl(N,Y,Z) 


N Ll 

N Li 

Nl 3 

n t 

Yl 3 

Yl 3 

Zl x 

Zl 7 

Zl 3 


PqHq 

pH 


(3.11) 


has non-singular leading order coefficients. This suggests that the St term in (3.96), which 
ensures the absence of singular points for St < 0, should be retained. Further, for St = 0 
we want to match the solution close to the stationary particle to a solution for y which 
is regular away from this point. Thus we want to retain those terms which ensure the 
absence of singular points away from particle at time St = 0, i.e. the l\- and /3-terms 
in (3.96) and the / 3 -term in (3.9a). The appropriate scaling is therefore 

/,=|«|i£i , l 2 = \St\*L 2 , l s = \St\H 3 , (3.12a, 6, c) 


n = |«|*JV , z = \6t\±Z , y=|«|-*y 


(3.12d, e, /) 


These scalings suggest that the separation process occurs in a relatively thin strip, n ~- 
|6t| 3 / 2 along a segment of the separation line x = T 0 (2,f) of length z ~~ |6/.| s. 

For the scaling (3.12), the solution for Y is most easily found by integration of (3.11) 
as in (3.9a), where L 2 and £3 are eliminated in favor of N and Z, which are constant along 


15 


the lines of integration, using (3.8). The result is: 


y _ ( t U iL + r L ° dL ) 

P.H, \J-oo </P(L\N,Z) J L, y/P(L-,tT,Z) ) 

where 

P(L-,N,Z) = -\nv{zln in L* + (3 n m Z 2 - 6 k{z\)L - 6 z\N) , 


(3.13a) 


(3.136) 


and Lo(N, Z) is the real root of the cubic P . This root is a unique and continuous function 
of N and Z since P is a monotonically decreasing function of L from (3.10). 

The choice of sign of the square-root in (3.13a), and the limits of integration are 
determined by the topology of the lines of constant N and Z. In physical space these lines 
are straight and pass vertically through the boundary layer; however, in Lagrangian space 
they are highly curved near the separation particle, as shown qualitatively in figure 2a. The 
lines can be divided into three segments corresponding to three asymptotic regions. The 
lower segments start at the wall and extend upward towards the vicinity of the separation 
particle. Because the Jacobian is nowhere singular along these segments, the t/-positions 
of the fluid particles remain finite on the boundary-layer scale, i.e. the scaled coordinate 
Y is small. Hence, these lower segments give rise to a layer of particles at the wall with 
a thickness comparable to that of the original boundary layer, this is shown schematically 
in figure 1. 

Along the central segments, the lines of constant N and Z pass through the vicinity 
of the separation particle. Here the y-position of the particles grows rapidly, and is given 
in scaled form by (3.13). Thus the central segments give rise to the intermediate, thicker, 
layer of particles shown in figure 1. The topology of the central segments in the Lagrangian 
domain, figure 2a, determines the choice of sign in (3.13a). From (3.8) and (3.9) it follows 
that on integrating upwards, Li increases from large negative values towards Lo(N,Z). 
Since Y is increasing along this part, the negative sign in (3.13a) applies. At position L 0 , 
the lines of constant N and Z turn around in the Lagrangian domain and Li again tends 
to — oo; along this second part the positive sign in (3.13a) applies. 

Along the third segments, the lines of constant N and Z proceed upwards toward the 
external flow. As in the lower segments, the Jacobian is no longer small here. Thus the 
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changes in y are finite on boundary-layer scale, and the third segments give rise to a layer 
of particles with a boundary layer scaJe thickness, atop the central region, as shown in 
figure 1. 

Hence, the separation structure is one in which the boundary layer divides into a 
central layer of physical thickness proportional to Re~*\8t\~* between two ‘sandwich’ 
layers of thickness proportional to Re~ *. 

The structure (3.13) is identical to the one obtained by Van Dommelen & Shen (1982) 
for two-dimensional separation, except that the coefficients now depend on the position Z 
along the describing line * 0 - A convenient way to illustrate the influence of the position 
Z is to scale out the coefficients using a procedure similar to Van Dommelen (1981): 


2 - 


i, =fci, =fc(z a + i)U; . + i)U; , (3.i4a,6) 


A r = a0\ N = a0\ (Z 2 + 1) * IV* , Y = — T = , Z = ^Z , (3.14c, d,e) 

7 /? 2 ’(Z’+l)i Pi 


where the tilde-variables scale out the Taylor series coefficients, the starred variables scale 
out Z, and 


<*=3*1111 1 7 =: (3 2 3 n 22**lll) 


\ P'H, 


PO s Hos ’ 


0o — 


z 3 n m 


=(^~y , 02=(~^ ± y 

\ z 3 n lu J \ tim / 


( z ! ti 2 2 tiiii ) j 


r > 0i 


In terms of the starred variables (3.13) reduces to 


r L ° dis r L ° 

J-oo Ji ;v* -3 T* - 1* 3 Jl : 


dL* 


where 


y/2N~ m -3 T* - I* 3 Jl\ y/2 N* - 3 T m - I* 3 


L- 0 (N-) = I(N') , 


and the function / is the inverse to the cubic N* = l/ 3 -f §/, i.e. 


(3.14d, e) 
(3.14/, g,h) 


(3.14t) 


(3.15a) 


1 1 
I(N*) = (N m + (1 + JV* 2 )>) 3 + (N* - (1 + N* 2 )^Y 


(3.156) 


The values of 0\/02, a 0 and J0 2 2 depend on the choice of the Eulerian coordinates (a:,z), 
but not on the definition of the Lagrangian coordinates. 
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An alternative expression for Y* can be found in terms of the incomplete elliptic 
integral of the first kind ): 




(3.16a) 


where 

/ \ i 1 Q T * 

A(AT*)=(3(I 0 * 2 + 1)) 4 , m(N*)= - + iX | , (3.166, c) 

^^.) = 2 arctan(^5) . (3.16d) 

Elliptic integrals are distorted identity functions, (in particular F(<p|0) = <p exactly), so 
that the arctan is responsible for the major variations in Y along the characteristics. 

Further terms in the asymptotic expansions (3.8) and (3.16) can be found in princi- 
ple. We note that the next term in the expression for Y does not involve a logarithmic 
correction, even though logarithmic second order terms do arise for the symmetric flows 
studied in the next section. 


We now turn to the physical interpretation of these results. The boundary-layer 
thickness is asymptotically determined by the position of the upper particle layer in figure 
1; letting L\ — ► — oo along the positive branch of (3.16a), we obtain the scaled boundary- 
layer thickness as 

r + ‘(A - *) ~ ■ < 3 - 17 ) 

The function Y + *(N*) gives the general shape of the boundary-layer thickness in a cross- 
section of constant z. For large values of N * the boundary-layer thickness decays toward 
zero much more slowly than suggested by the sketch in figure 1. Nevertheless, at the outer 
edges of the thin separation region, the solution still matches with finite values of y; for 
from (3.12), (3.14), and (3.16) 


4al F (l 

3« 2s y V 2 


1 \/3 

2 ± 4 


In |* 


for |6t| J -C |n| -C 1 


(3.18) 


To show the dependence of the boundary-layer thickness on the coordinate z, contours 
of constant Y + in the iV,Z-plane are plotted in figure 2b. Note that the coordinate N 
is measured from the oblique, curved, separation line. Actual lines of constant boundary- 
layer thickness might, for example, appear as sketched in figure 2c, which has been drawn 


18 



by taking |$<| = 0.06 and unit values for various coefficients in (3.7) and (3.14). Asymptot- 
ically, the boundary-layer thickness has the shape of a crescent shaped ridge. The crescent 
shape is long and thin, i.e. quasi-two-dimensional, because from (3.12) the n length scale 
is asymptotically shorter than the 2 length scale (note that for three-dimensional steady 
separation Smith (1978) has proposed a quasi-two-dimensional structure). In an Eulerian 
numerical calculation, the development of such a crescent-shaped ridge may be a possible 
diagnostic indicating the presence of a singularity. 

Evidence of this type of singularity is provided by Ragab’s (1986) calculations for 
impulsively started flow past a 4:1 prolate spheroid inclined at a 30° angle of attack. His 
results strongly suggest that the displacement thickness becomes unbounded away from 
the symmetry line. However, it is not possible to deduce the shape of the singularity from 
the results presented. 

A point of interest is the decay of the boundary-layer thickness along the describing 
line for large Z. From (3.12), (3.14), and (3.16), 



for |6t| * <C \z\ < 1 


(3.19) 


Hence for increasing z, the separation structure expands in n-direction, while the thickness 
of the boundary-layer decreases. 

The particle propagation velocity n which gives rise to the accumulation of particles 
at the separation line is, according to (3.8a), given to leading order by 


n ~ |£t| 2 n\Li 


(3.20) 


To describe this in the more familiar Eulerian coordinates, the transcendental relationship 
(3.16) must be inverted to the form 

L\ =L\{N*,Y '*) . (3.21) 


The inversion has been performed numerically, and in figure 2d we present contours of L\ 
in the (A 7 *,K*) plane. From (3.14), 


n — 


■i«i s + 1)‘£; + + 1}i ) ' 


( 3 . 22 ) 
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It follows that the lines of constant L\ shown in figure 2d describe the shape of the lines 
of constant n in cross-sections of constant ~z through the separation structure. They also 
give the asymptotic shape of the lines of constant velocity components x and z and density 
p in these cross-sections, since 

(x,Z,p) = (x s ,Z 3 ,p s ) + | 6 <| ’ ^(x 1 ,Z 1 ,p 1 )P 2 (Z 2 + + (x3,Z3,p 3 )-j^Z^Z^ + . . . 

(3.23) 

We note that the topology of figure 2d for ~ 0 seems quite close to the computed lines 
of constant velocity presented by Van Dommelen (1981) for finite |6<|, and thus lines of 
constant velocity might be a useful indication of an incipient unsteady separation. 

The next point of interest is the shape of the velocity profiles. According to (3.23), 
in Eulerian space the velocity profiles must develop a large flat region of nearly constant 
velocity as separation is approached. However, accepting the numerical results of Van 
Dommelen (1981), this flat region is only evident extremely close to the singularity, so that 
resolution problems or finite Reynolds number effects tend to obscure the phenomenon. 
From (3.23) and figure 2d, the velocity profiles near an incipient three-dimensional separa- 
tion must have a local maximum or minimum in velocity. However, this is not necessarily a 
precise indication of incipient separation. For example, in the case of the circular cylinder, 
a minimum in the velocity profiles develops relatively quickly, after | diameter motion, 
yet separation occurs much later, after | diameter motion. Figure 2e shows the shape of 
the velocity profiles near the interior extrema. The shapes of the velocity profiles in the 
sandwich layers at the edges of figure 2e cannot be found from asymptotic analysis since 
they depend on the precise details of the earlier evolution (cf. the remarks below (3.24) 
and part 2). 

A more significant sign of the start of separation might be a transverse expansion 
of the lines of constant vorticity near the velocity minimum/maximum; since the above 
analysis is inviscid to leading order, the vorticity lines closely follow the motion of the 
boundary-layer particles. In the boundary-layer approximation, the vorticity is the y- 
derivative of the velocity distribution. The corresponding asymptotic topology of contours 
of dL*/dY * is shown in figure 2f. This topology seems close to the computed vorticity 
lines presented by Van Dommelen (1981) for a time near separation. 

The asymptotic structures of the upper and lower vorticity layers are similar to the 
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two-dimensional case (Van Dommelen 1981). Expressed in terms of Eulerian coordinates, 
they take the form of regular Taylor expansions: 

(*,!/, z,p)= X ( u Htnr(y),v~ nr (y),w~ nr (y),p~ nr (y)^ , (3.24a) 

mnr>0 

and 

^ ^ ?ft, (4nr(j).*L(y))®mr(j).Pinr(s)) . (3.246) 

mnr>0 

respectively, where the sums run over the non-negative integers, and the Prandtl transfor- 
mation y = y — y* (x,z,6t), describes the motion of the upper layer. 

Substituting (3.24) into the boundary-layer equations, we find that the u^ nr . tu^ nr , 
Pmnr> > 0, v > 1) and the t^ nr , (m,n,r > 0) are determined in terms of the 

(^mnoi^mno^mno)? but that these latter functions axe indeterminate due to the depen- 
dence of the solution on earlier times. The (u^ n0 ,wtino ,P™ n0 ) must, however, satisfy the 
boundary conditions (2.5a) at the wall, and match both at the outer edge of the boundary 
layer (see (2.56)), and with the central inviscid low-vorticity region. At fixed N and Z, 
the latter matching conditions yield from inverting (3.16) and using (3.23), 

^^ 41 

( u ooo 7^000 > Pooo ) * iPs ) “ (®i ? z \ iPi )~ 2 35 2/ ¥ +00 , (3.25a) 

7 V 

_ 4 1 

(«ooo i™ooo »Pooo ) (*3,Zs,Ps) ~ (xuz u p t ) — — as y -> -oo . (3.256) 

7 y 

Asymptotic matching conditions can also be derived as |JV|,|Z| — ► oo, as Van Dommelen 
(1981) has done for two-dimensional flows. 

A final point of interest is the ‘accessibility’ of the region of flow beyond the time 
of initial separation. In a steady Eulerian computation, Cebeci, Khattab & Stewartson 
(1981) took the accessible region to be the domain where a boundary-layer solution can be 
found, (whether it is still an asymptotically correct solution of the Navier-Stokes equations 
in the presence of interaction or not). In the Lagrangian case, some care is needed, because 
the singular continuity equation is integrated separately. Numerical experiments such as 
the one in part 2 do in fact suggest that the non-singular momentum equations can be 
integrated past the separation singularity without apparent difficulty. When that is done, 
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the vertical lines through the boundary layer appear in Lagrangian space as shown in figure 
2g rather than figure 2a. For the shaded particles, y is indeterminate; these particles may 
be thought of as having disappeared at infinite y. Yet the continuity equation can still 
be integrated along all lines of constant n and 2 which start at the wall. A singularity 
develops only on the line passing through the saddle point in figure 2g, which for 0 < St <C 1 
corresponds to a singular line segment 



(3.26) 


However, the solution so obtained must be considered meaningless at least for all particles 
which have at some previous time passed through the singular curve. For that reason, we 
define the region of inaccessibility as those stations (*, 2 ) which contain particles which 
have at any time been on the singular curve. Initially, the region of inaccessibility will 
primarily expand in the 2 -direction through the scaling (3.12e). In the n-directior it will 
expand by means of the motion of the describing line (3.7) and additionally through the 
motion of the particles which propagate downstream away from the singular curve. Thus, 
the region of inaccessibility extends over a finite surface area, rather than just the curve 
(3.26), in agreement with the steady Eulerian definition of Cebeci et al. (1981). 

Naturally, the singularity structure derived here will not remain asymptotically correct 
arbitrarily close to t = t s , because the normal velocity above the central inviscid region 
becomes infinite at t — t s . From a study of the Navier-Stokes equations it is found that the 
singularity is smoothed out when a ‘triple-deck’ interaction comes into operation for St = 
0(Re~" ), at which point the scaled boundary-layer thickness is 0(Re ”). Because the 
singularity is quasi-two-dimensional, the scalings and governing equations axe essentially 
those derived by Elliott et al. (1983) for two-dimensional flows, but with the addition of a 
passive 2 -momentum equation. In the central interaction problem, the coordinate 2 , which 
has an interaction length scale 0(J?e _TT ), only appears as a parameter. However, it is not 
clear whether the singularity will be completely removed by the interaction (Smith 1987). 


4. Three-dimensional symmetric separation 

In the previous section a separation singularity structure was derived assuming that 
the flow was arbitrary, an assumption that might be appropriate for flow past an asymmet- 
ric body. However, in the case of a spheroid at relatively small angles of attack it is likely 


22 



that separation first occurs on one of the symmetry lines; indeed numerical calculations 
confined to the symmetry line have been performed on this basis (Wang & Fan (1982) 
Cebeci, Stewartson & Schimke (1984)). In sub-section (a) below we derive the form of the 
singularity appropriate for separating flows where the separation line crosses a symmetry 
line normally. 

However, this is not the only type of symmetric separation of interest. When a sphere 
is impulsively rotated about a diameter, centripetal effects generate a boundary-layer flow 
towards the equator. After a finite time an equatorial singularity develops as a result of 
a boundary-layer collision. The structure of this singularity on the symmetry line has 
been determined by Banks &: Zaturska (1979), and Simpson & Stewartson (1982a). In 
this case the separation line coincides with the symmetry line. Similar singularities occur 
after a finite time at the apex of a horizontal circular cylinder which is impulsively heated 
(Simpson Sc Stewartson 1982b), at the inner bend of a uniformly curved pipe through which 
flow is impulsively started (Lam 1988), and at the stagnation points on a two-dimensional 
cylinder in oscillating flow as a result of steady streaming effects (Vasantha Sc Riley 1988). 

A more general form of the singularity generated by two symmetric colliding boundary 
layers on a smooth wall would first develop at a point rather than along the entire symmetry 
line. For example, such a singularity might develop on the equator of an ellipsoid which 
is rotated about one of its principal axes, or in starting flow through a curved pipe with 
non-uniform curvature, or at the apex of a heated ellipsoid. In sub-section (b) the three- 
dimensional structure of such a singularity is derived. The results on the symmetry line 
agree with previous authors, but the simplicity of the Lagrangian approach allows us to 
determine additionally the singularity structure off this line. The latter is a necessary 
preliminary in order to formulate subsequent asymptotic stages in the separation process. 

Another class of separation singularities are rotationally symmetric about the sepa- 
ration point, so that the separation line degenerates to a point. For example, singularities 
develop after a finite time on the axis of a spinning disc or sphere whose direction of 
rotation is impulsively reversed (Bodonyi & Stewartson 1977, Banks Sc Zaturska 1981, 
Stewartson, Simpson Sc Bodonyi 1982, Van Dommelen 1987), and at the apex of a sphere 
which is impulsively heated (Brown Sc Simpson 1982, Awang Riley 1983). The structures 
of these singularities, which differ due to the presence and absence of swirl, are derived in 
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sub-sections (c) and (d) respectively. The results on the axis agree with those of previous 
authors, while the singularity structures off the axis are new. 

(a) Lateral symmetry 

When the boundary-layer flow is symmetrical about a line along the surface of the 
body, the describing line of separation must either cross the symmetry line normally or 
coincide with it. In this sub-section we will address the case of normal crossing, leaving 
the second possibility to the next sub-section. 

For consistency with section 3, we identify the compressed coordinate n with * and 
take the (, T 7 -plane as the symmetry plane so that a: is an even function of £ and z an odd 
function. Then the analysis is a simpler version of the one in the previous section. The 
only transformation of the Lagrangian coordinate system needed is a rotation around the 
£-axis to eliminate the x J2 derivative. Also, the discussion concerning which derivatives 
must be zero if t 3 is the first separation time (see (3.4) and following) can be restricted to 
the symmetry plane to show that the second order derivative which is forced to be zero 
must lie within the symmetry plane. 

Hence the structure of the separation process remains basically unchanged, although 
the describing line of separation simplifies, and is now symmetric about the symmetry line 
z = 0 (cf. (3.7)): 

n = x - x(£ 3 ,t) - ^jz 2 . (4.1) 

A degenerate case is two-dimensional flow, where x is totally independent of £, and the 
separation line becomes a straight generator in the z-direction. In addition, the coefficient 
0i vanishes, which suppresses the decay of the boundary-layer thickness with z. The 
resulting structure is described in detail by Van Dommelen (1981). 

Thus lateral symmetry, or more strongly two-dimensionality, does not fundamentally 
alter the separation process. This conclusion is consistent with the symmetry line calcula- 
tions of Cebeci, Stewartson & Schimke (1984). 

(b) Symmetric boundary-layer collision 

When the describing line coincides with the symmetry line, significant changes in 
structure are unavoidable, since the flow is symmetric while the separation structure illus- 
trated in figure 2 is asymmetrical. 
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We identify the compressed coordinate n again with ar, but now we assume that the 
T 7 ,£-plane is the symmetry plane, so that x is an odd function of £ while z is an even 
function. A singularity occurs when x^ first vanishes at the symmetry plane, since the 
derivatives x tV and x t ( are zero by symmetry. Since the first occurrence of a zero value 
must occur where x^ is a minimum, the second order derivatives X( v and x ^ must vanish, 
while the other second order derivatives are zero by symmetry. 

The fact that all the second order derivatives are zero invalidates the scalings for 
t] and y made in the previous section (e.g. (3.12), (3.14)), hence a separate analysis 

with significant modifications is needed. Proceeding along similar lines as in the previous 
section, a local Lagrangian coordinate system ki,k 2 ,lc 3 is introduced with origin at the 
separation particle, but with the same orientation as the original axis system. A rotation 
of this coordinate system around the fcj-axis, 


fci = ( , k 2 = 


23 h - z 2 k 3 


v 7 z l + 4 ’ 3 V z l + 4 

X ~ x , z = z-z(( a ,t) , 

can be made to eliminate the z 2 -derivative. The shearing transformation 


fca = 


z 2 k 2 + Z3 &3 


h — £ > l 2 = k 2 + -zr £3 > h — ^3 j 

®122 


® = ® , z = z-z(t,,t) , 


eliminates the *123 derivative, resulting in the Taylor series expansions 


(4.2a, 6, c) 


(4. 2d, e) 

(4.3a, b, c) 
(4.3d, e) 


* ~ + 2 * 122^2 + 2* 133 ^3 + •••+ StX\l 2 + ... , (4.4a) 


Z ~ Z3I3 + . . . 


(4.46) 


The expressions for the characteristics of the Jacobian equation for y become 


dh 

dy 


pH 

PoH 0 


h { — 2 3®122 ^2 + • • •} 


( 4 . 5 a) 


dl 2 
dy 


TT 

{ 2 3(|® llll? + §®122^2 + 2*1 33 ^3 + ^®l) + •• •} 


(4.56) 
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In order to avoid singularities for Si < 0, the quadratic in (4.5b) must be elliptic and of 
opposite sign to SFj. Sinre x ^ is initially positive, rf. (2.1), it follows from (4.1a) and (4.5b) 
that at a first zero 


*in > 0 i *122 > 0 , *133 >0 , *i < 0 . 


(4.6a, 6, c, d ) 


The topology of the characteristics (4.5), shown in figure 3a, can be compared to the 
asymmetric case figure 2a, where the separation characteristic develops a cusp at St = 0. 
In this case, the separation characteristic is constrained by symmetry to remain straight. 

Appropriate local scalings near separation can be found using arguments similar to 
those of the previous section: 


h = m-'fa'L, = \6t\*fc(Z 2 + l)*Li , 

h = = |«|*AA(Z 2 + 1)U; , 

* = \6t\’a/3l X = \St\ 3 *a0l{Z 2 + 1)*A~ , 

Y 


Y* 


r . t = , 


|6<|*7/?2 \St\\^ 2 {Z 2 + l)i ’ ' ' /?i 

Psffs 


a = 5*m i 7 = (5*3*122*111) 


po,H< 


Oj 


/?o 


*3*111 R — ( g !33 \ 2 a _ ( ^ X 1 \ 

' (^3*122 *11 1 ) * ’ ^ ~ \ *3 *1 1 1 / ’ ^ “ V *11 J 


The continuity integral becomes 


f L 0 

1 


dL* 


f L ° d£_ 

Jl ; JL*(2X* - 31* - 


i/i*(2A r * — 3L* - IF) J l\ s/TF{ 2 


I* 3 ) 


where 


LW) = /(A*) . 

This can be written as an elliptic integral similar to (3.16), 




(4.7a) 

(4.76) 

(4.7c) 

(4.7<f,e) 
(4-7/, 9) 

(4.7 h,i,j) 


(4.8 a) 
(4.86) 


(4.9a) 
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where 


a(x-) = (3<l; : -h)(i; j +3))‘ 


m , Y .v i u? + 6 

m(JC > - 2 “ ~Ta5“~ 


vK£;,JT*) = 2arctan ^ ( L f + 3 ) (f| “ J ) ) ' 


(4.96, c) 
(4.9(f) 


Note that instead of using TJ as the independent variable, there is an advantage in using 
L 2 , as given implicitly by the relation 


£;(£;, j n = (if + 1 )\hx-/(l? + 1)>> , 


(4.9e) 


since at the symmetry line the solution is regular in terms of L 2 : 

2 


r*(I 2 *,0)~-^(^+arctanI*) . 


(4.10) 


Contours of the boundary-layer thickness Y + in the X,Z plane are shown in figure 
3b. The asymptotic relations for large |A r | and \Z\, corresponding to (3.18) and (3.19), are 


4a 3 

1 Vs\ 1 

3i2^ 7 V2 

2 4 / |*| 3 

„ 1 w 

* \ 1 

1lh Y l 

a/3j 3 |z 3 / z 


|*| < 1 


(4.11) 


(4.12) 


The velocity components and density in the neighbourhood of the stationary point 
are given by 

i~-|«|)§a(J|(Z J +l)U; , (4.13a) 

(i^)~(i„ P .) + |<(|^(i 1 ,? 2 )AA(Z J + l)Si;+(j s ,?,)^-Z^ . (4.136) 

Hence L\ can be interpreted as the velocity component, towards the symmetry plane. The 
scaled velocity profile, — LJ, is illustrated in figure 3c at a number of A'* stations, while 
contours of L\, and the corresponding vorticity component, dL^/dY*, are illustrated in 
figures 3d and 3e respectively. In cross-sections of constant z, the variations in velocity 
parallel to the symmetry plane are proportional to L 2 . L 2 velocity-profiles are given in 
figure 3f, while figures 3g and 3h illustrate contours of L 2 and the vorticity dL 2 /dY*. 
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Close to the wall, i.e. as Y — > 0, 

i~-j7 2 #*V 2 , (*,/>) ~(w#)“ (i 2 ,p 2 )~ — ~ , (4.14a, b) 

4 V3 7 y 

which match to a regular vorticity layer of similar form to (3.24). Similarly a match can be 

achieved with a separating layer governed by a Prandtl transformation above the central 

inviscid region. 

Figure 3i shows the characteristics, i.e. lines of constant x and z, for 6t > 0. Inte- 
gration of the continuity equation yields a singularity over a segment — 1 < Z < 1 of the 
symmetry line X = 0. Since neither the singularity nor any particles on the symmetry line 
leave the symmetry line, the region of inaccessibility remains restricted to the symmetry 
line. 

A special case occurs for separation at the intersection of two symmetry lines, such 
as at the apex of an ellipsoid. In that case, in addition to the symmetry in £, x is an even 
function of £ and z an odd one, and the transformations of the Lagrangian coordinate 
system (4.2) and (4.3) become trivial. No changes in the leading order singularity structure 
occur, since it was already symmetric in z-direction, even though this condition was not 
imposed. However, the velocity parallel to the symmetry plane must be antisymmetric, 
and the density symmetric (cf. (4.13b): 

— Z ’ P~P, + MWo/M£ 2 +l)U; • (4.15a, 6) 

Pi *3 

In the case of two-dimensionality, where a; is independent of £, the coefficient 3\ van- 
ishes as in the previous subsection, suppressing the decay of the boundary-layer thickness 
with z. The flow on the symmetry line can then be written as a one-dimensional problem, 
and was studied from an Eulerian standpoint by Banks and Zaturska (1979) and Simpson 
fc Stewart son (1982a,b). In part 2, Van Dommelen (1989) uses this flow to verify the 
Lagrangian analysis numerically to high accuracy. Favourable numerical comparisons with 
the singularity structure away from the symmetry line have been obtained by Lam (1988) 
for starting flow through a circular pipe. 

The existence of this singularity has also been reported by Stern Paldor (1983), 
Russell Landahl (1984) and Stuart (1987) while studying inviscid models for the growth 
of large amplitude disturbances in boundary layers. In fact, because unsteady separation is 
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primarily inviscid in its final stages, an alternative approach to that above would be to solve 
the inviscid version of (2.3) exactly, and then to examine the possible singularities of the 
solutions (see also Van Dommelen (1981) for the two-dimensional singularity). Note that 
although Stuart’s (1987) exact, inviscid, symmetry line solutions do not include a parallel 
flow in the z-direction, our results are in agreement since the details of the singularity 
structure are independent of z(£ 4 ,f). 

As in section 3 the above singular solution will not remain valid for sufficiently small 
|6t| because previously neglected pressure gradients will become important (cf. the inter- 
active problem for the two-dimensional singularity formulated by Elliott et al. (1983)). 
Further, because the velocity towards the separation line is much smaller in the upper and 
lower vorticity layers than in the central layer, it is in the vorticity layers that the effect of 
the pressure gradient will be felt first. However, it is the central layer which is responsible 
for the growth in boundary-layer thickness; thus it appears that the first asymptotic rescal- 
ing does not lead to an ‘interactive’ effect to smooth out the above singularity. Instead, the 
singularity continues to be driven by the flow in the central layer, while significant changes 
occur in the upper and lower layers. Similar arguments seem to hold for the singularities 
in (c) and (d) below. 

(c) Axisymmetric boundary-layer flow with swirl 

In axi-symmetric flow, the flow geometry does not depend on £ and £ individually, 
but only on the Lagrangian distance, 

0 = V? + < 2 , (4.16) 

from the axis £ = ( = 0. The displacement of rings of particles g = rj = t = constant from 
their original position must remain restricted to a change in physical distance, 

r = , (4-17) 

from the axis, a rotation around the axis, and a shift in vertical position. Hence according 
to the theory of orthogonal matrices, the solution must be of the form 

* = c(p 2 ,T7,<)£ + j(p 2 ,»M)C * (4.18a) 

2 = s(q 2 ,Vi + c(0 2 ,V, 0C , (4.186) 
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(4.19a) 


where, because of the assumption of regular * and z, 


c ~ ® i{ (0,T7,0,<) + !*,«{(0,T7,0 ,t)g 2 + ... , 

s ~ ® jC (0,T/,0,t) + |x )CC< (0,7/,0,t)e 2 + ... . 

In terms of c and s the physical distance from the axis is given by 


r 2 = (c 2 + s 2 )g 2 


The Jacobian J in (2.2) can be written in terms of g and r as 


3 = ( r2 ),e , y.t - • 


(4.21) 


Thus separation occurs at a stationary point for r 2 (g 2 and from (4.20) and (4.21) it 


occurs on the axis when 


c(0,T] 3 ,t a ) = s(0,77,,t a ) = 0 


(4.22a, b) 


A rotation of the Lagrangian coordinate system to diagonalize the second order deriva- 
tives of x is not advantageous here, since the axial symmetry would be lost. Instead we 
rotate the coordinate system around the symmetry axis, 

f *23^ - C r 7 *12 £ T ®23 C L \ 

kl = / 2 T ’ = V - V3 , *3 = -r = f= , (4.23a, b, c) 


*12 + *23 


*12 + *23 


to eliminate the ®i 2-derivative, followed by the shearing transformation 

h = ki , h = k 2 + jr^-(fci + kj) + St^~ , l 3 =k 3 , (4.24a, b, c) 

0®23 *23 

to eliminate 2:333 and *3. 

The characteristics of the Jacobian are lines of constant distance r from the axis. If 
t s is the first time that a singularity forms then ®m*i must be negative, or for a suitable 
choice of the positive /j -direction 


*111 > 0 T j < 0 


(4.25a, b) 


The characteristic lines of constant r in the g,l 2 -plane appear as sketched in figure 4a. 


30 



Appropriate local scalings are 


Y* 


e = \St\>0P' , (, = \6t\M 2 L; , r = , y = , (4.26a, b, c,d) 


a = 3*111 , 7 = 3*111 (® 33 ) * rj’ , A - J 11 ! , 0 - (-=— ^) , ( 4 . 2 6 e,/, 0 , 

POj-^Oj 2( 5 C 23 ) 2 V ®111 / 

leading to a continuity integral 

2dP* 


h) 


Y m ~ f 

Jo 


p ° 2dP* 


r 


^/4ii* 2 — (3p* + p* 3 ) 2 y P . v 4 ^* 2 - ( 3P * + P * 3 ) 2 ’ 


where 


p;(/r) = /(ir) . 

This can be written as the elliptic integral 

}'•(!;, R-)= |F(j|m)±if(»>|,n) , 

where 


*>(£;.*•) = 2 arctan (p; 3 + 3) ’ - l)) , 

and P* 2 is related to Lj and R* through the solution of the cubic equation 


4 R* 2 = 9i; 2 P* 2 + (P*^ 3) 2 p* 2 


(4.27a) 


(4.276) 


(4.28a) 


A(R-) = (3(P; , +l)(Po‘ 2 +3) 3 ) i , m (R-)=l- ^” 4+ ^p 2+18 , (4.286, c) 


(4.28d) 


(4.28e) 


On the axis, (4.28a) simplifies to 


r*(I 2 *,0) = 1(1 + arctan , 


while for large R * , the bound ary- layer thickness asymptotes to 


, 4oa / 7T 1 v/3\ 1 r , 

y ~ 3$2iy F \2 2 ~ T )""* 


2 

r 3 


(4.29) 


(4.30) 
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The velocity components in the radial and azimuthal directions, and the density, are 


r ~ +xw;) , re ~ - x u;) , * = ^ , (4.3ia,t,c) 

X 1 






P = P. + \6t\(p 2 /3 o 0 2 L; + \p u 0* P* 2 - p) , 


where c = sgn(«23)> and 


. _ (P* 2 + 3)P* 2 3P* 2 Z-2 




(4.31rf) 


(4.32a, b) 


W* — 

2R* ’ 0 2R* ' 

are the symmetric and anti-symmetric velocity profiles shown in Figure 4b. 

Both the radial and the circumferential velocity profiles depend non-trivially on the 
parameter \ • However the magnitude of the velocity, 


q — \/u 2 + w 2 = |£f| * §a /? 3 y/l + x 2 P* , 


(4.33) 


does not; contours of q are illust rated in figure 4c. The vorticity components normal to 
the velocity and parallel to it are proportional to 


u;u; Y . + w;w; Y . 3 , 

~' 5 J ’ 


n p = 


w^y. - u;w 0 \. 3 p.,p.2 


p*(p* 2 + 1) . 


(4.34a, 6) 
(4.35a, b) 


VUp+W 

Contours of these quantities are plotted in figures 4d and 4e respectively. 

A match with the sandwich layer adjacent to the wall is again possible, since as 1' — ► 0 


2 l 

r ~ \lfP 2 xry . r 0 ~ -\nP 2 <;ry , p~p,-p 2 - — - 

3 7 y 


(4.36a, b, c) 


Similarly a match can be achieved with the upper separating layer. 

Figure 4f shows the characteristics for St > 0. The singular line is the physically 
expanding circle R* = 1, but the region of inaccessibility is larger due to particles with 
^ 0 which move radially outward from the singular line at a greater rate. 

On the axis itself, the continuity integral is particularly simple: 


_ [ v ppHpdn 

y Jo pH (*?( + *? c 


) ' 


(4.37) 
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When this integral is expanded to second order, a logarithmic correction to the y-position 
is obtained. This and other terms were initially overlooked in Eulerian analyses of the flow 
on the axis (Bodonyi & Stewartson (1977), Banks & Zaturska (1981)), but in a Lagrangian 
approach the logarithmic term follows naturally from the hypothesis that the solution for 
the motion parallel to the boundary is regular. (A similar logarithmic term arises in the 
symmetric case (b) above, cf. part 2). In fact, from this hypothesis alone, the complete 
singularity structure presented by Stewartson, Simpson & Bodonyi (1982) can be recovered 
by means of a simple integration of (4.37). 

(d) Axisymmetric boundary-layer collision without swirl 

Finally, we consider the case of axially symmetric flow when there is no rotation of 
the flow about the axis. In the absence of such rotation (4.18) simplifies to 


* = * = i’MK 


(4.38) 


No transformations of the Lagrangian coordinate system are needed in this case. It follows 
that if a singularity first appears on the axis c must vanish. The contours of constant r 
are then identical to those for a symmetric collision (figure 3a), while the Taylor series 
coefficients satisfy conditions (4.6a,b,d). 

In a similar way to before suitable local scalings are 

0 = \6t\i(3P* , k = \6t\ l ’/3 0 /3L ; , r = \ 6 t\l a /l 3 R* , y = • - , (4.39a, 6, c, d) 

M 5 iP 


, 7 = & = (r L y ’P-l-T 1 -)’ ’ 

pOsttos \Xl22 / \ * 111 / 

(4.39e,/,<7,/i) 

leading to a continuity integral 

r p ° p**dp* ( r p '° p*»dp* 

y 0 R*y/ 2 R* -3P* -P* 3 ± /p* R*y/ 2 R* - 3P* - P* 3 ’ ' ^ 

where 

P ;(R-) = I(R-) . (4.408) 


This solution can be ‘reduced’ to the form 

1 / 4 - 4n 


y*(l;,r*) 


p; + 3 


. . 7r* . . 2 . 7T . . 

B(n; -\m) - -F(- |m) 
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± i^(^ n(n; * |m) -X FMm) 


2 fK sin(^) \ 

Po \2A y/l — msin 2 y?/ 

2 4 f nP* 0 sin(2y?) 

+ — axctan — — 7 . ■ = = 


2A wT- 


msin ip 


(4.41a) 


where II(n;y?|m) is the incomplete elliptic integral of the third kind, and 

1 1 3P* 2 + 6 

A (JT) = (3(P 0 * 2 + 1)(P; 2 + 3)) 5 , m(R') = - - , n (JT ) = 1 


x\/( p o 2 + 3 ) (f - 1 


<p(LZ,R*) = 2arctan I — W (PJ 


p*(i;,/r) = (£; 2 + i)>i(r*/(l; 2 + i)>) 


3 p ; 2 + 12 

2A 2 + 6 ’ 
(4.416, c, d) 

(4.41e) 


(4.41/) 


On the axis (4.41a) simplifies to 


Y*{L;,Q) = 4- axctan L\ + 


L? 4- 1 


while for large R* the boundary-layer thickness asymptotes to 


tiH 1 - 


^3. 7T 1 _ \/3 

T'2 2~ 4 




The velocity and density are given to leading order by 


(4.42) 


(4.43) 


f ~ -|6t|’§a/? 3 P* , p~ Ps+M'ipiML* 


(4.44a, 6) 


Sample velocity profiles are illustrated in figure 5a, while contours of P * and the vorticity 
dP * /dY* are given in figures 5b and 5c respectively. Again, a match is possible with the 
wall layer, since as Y — * 0, 


35-22 ,-2 3 A 1 

r ~ T 0'y>ry> , p ~ p, + P 2 ~T ~T ~T 

23 3 « 7 * y ? 


(4.45a, 6) 
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5. Discussion 


In this paper we have shown that the description of attached flow past a body using 
the classical boundary-layer equations can break down after a finite time due the formation 
of a local singularity. In a Lagrangian description the class of singularities is characterized 
by a singular continuity equation, but a regular momentum equation. The evidence shows 
that such singularities axe both mathematically consistent and physically relevant (e g. Van 
Dommelen k Shen 1980a, Van Dommelen 1987, 1989, Lam 1988, Bouard k Coutanceau 
1980). The precise structure of the singularity depends on the symmetry of the flow, 
and some of the simpler structures have previously been partially or totally described in 
Eulerian coordinates by other authors. The purpose of this paper is to provide a unified 
theory to facilitate the identification of singularities of the Lagrangian type when they do 
occur. This seems especially relevant for the difficult problem of the asymmetric singularity, 
where the singular behaviour would have to be deduced from a three-dimensional unsteady 
computation. 

These singularities occur when a fluid particle becomes compressed in one direction 
parallel to the boundary. Conservation of mass then implies that the fluid above this fluid 
particle is forced out of the boundary layer in the form of a detached vorticity layer. A 
common feature of all the singularities is that the typical length scale in the direction 
of compression is 0(|6t|*). However, the the strength of the singularity increases with 
the symmetry of the flow; the the boundary layer thickness varies from 0(|6t|~i) for the 
asymmetric symmetry to C>(|£t| - *) for the axisymmetric singularity without swirl. 

Because the singularities take the form of a vertical ejection of fluid from the boundary 
layer, we believe that they indicate the onset of separation as hypothesized by Sears k 
Telionis (1975). While the present singularity structures do at least seem to describe the 
initial genesis of the separating shear layer, within an asymptotically short time interactive 
effects which axe neglected in the classical boundary-layer formulation must be included 
(e.g. Elliott et al., 1983, Henkes k Veldman 1987). At that stage a new asymptotic scaling 
must be substituted into the Navier-Stokes equations in order to recover the correct large- 
Reynolds-number limit. Knowledge of the precise asymptotic structure of the singularities 
is necessary to identify this new scaling, and one of the contributions of this work has been 
to identify the full structure of a number of symmetric singularities. 
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At first sight the symmetric singularities may appear less likely to occur in problems of 

practical importance. However, they have previously arisen in inviscid models of ‘transition 
to turbulence’ in regions where symmetric counter-rotating longitudinal vortices are forcing 
the convergence of fluid particles (Stern et al. (1983), Russell et al. (1984), Stuart (1987)). 
A three-dimensional extension of the work by Smith & Burggraf (1985), may lead to 
an asymptotic description of transition which accounts for viscosity, where the turbulent 
bursts are associated with local regions of classical boundary-layer separation (symmetric 
or otherwise). 

This work was presented in part at the IUTAM “Fluid Mechanics in the Spirit of G.I. 
Taylor” Conference, April 1986, Cambridge. The authors acknowledge financial support 
from the NASA “Materials Processing in Space” program, the SERC, the AFOSR, and 
ICOMP, NASA Lewis. 
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toward the 
external flow 


toward the wall 

(A) LAGRANGIAN TOPOLOGY OF VERTICAL LINES THROUGH THE BOUND- 
ARY LAYER NEAR THE SEPARATION PARTICLE. 



(B) CONTOURS OF THE SCALED 
BOUNDARY-LAYER THICKNESS 
Y + = 5 1/2, 5, A 1/2, 

. . . , 2 1/2 IN SCALED, 
OBLIQUE COORDINATES. 



(C) POSSIBLE ACTUAL APPEARANCE OF CONTOURS OF BOUNDARY-LAYER 
THICKNESS (SCHEMATIC). 



(D) CONTOURS OF THE SCALED VELOCITY L-| = 0, ±1, ±2... IN 
SCALED COORDINATES. 

FIGURE 2. - STRUCTURE OF ASYMMETRIC THREE-DIMENSIONAL SEPARATION. 













(A) LAGRANGIAN TOPOLOGY OF CB) VELOCITY PROFILES OF THE TWO COMPONENTS -Uq AND wj 

PHYSICALLY VERTICLE LINES. 



IV 


(C) CONTOURS OF THE SCALED ABSOLUTE VELOCITY P* = 0, 1/2, 1, 

1 1 / 2 ,... 

FIGURE A. - STRUCTURE OF AXI-SYMMETRIC SEPARATION WITH SWIRL. 
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(F) LAGRANGIAN TOPOLOGY OF PHYSICALLY 
VERTICAL LINES BEYOND THE FIRST 
SINGULARITY. 


FIGURE A. - CONCLUDED. 
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